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1. INTRODUCTION 

In Morocco, the production of renewable energies is growing faster than that of fossil fuels. This 
dynamic should continue in the future thanks to the significant national and international investments which 
will allow the country to gradually free itself from its energy dependence both in terms of its imports of 
hydrocarbons and electricity imported from abroad. As part of Moroccan energy transition strategy, wind 
energy is recommended in many scenarios to achieve the ambitious objectives set up, compared, 
economically and technologically, to other sectors of electricity production of renewable origin [1], [2]. 

In the literature, to control the active and reactive powers of the double fed induction generator 
(DFIG), the authors [3], [4] have respectively proposed the use of direct power control (DPC) and oriented 
stator flux vector control. Active and reactive power control is obtained with a controller called rotor current 
regulator or by modified direct power control (MDPC) [5]. However, the wind conversion system (WCS) is a 
complex system with non-linearities, multiple variables, and system parameter disturbances caused by the 
external environment. These problems should be taken into account to meet the challenge to obtain an 
accurate mathematical model [6]. Therefore, designing more detailed and effective control methods to solve 
these problems of different disturbances affecting the system is of great importance for the reliable operation 
of the WCS and also for extracting the maximum power from the wind turbine. Currently, many researches 
have been presented [3], [4] using proportional integral (PI) regulators and backstepping control [6] to 


Journal homepage: http://beei.org 


Bulletin of Electr Eng & Inf ISSN: 2302-9285 D 1211 


control rotor currents. But the problem is that its last two regulators are sensitive to the variation of the 
internal parameters of DFIG. Some authors have looked for other power control options for DFIG using 
model predictive control [7], sliding mode control [8], [9] or feedback linearization control [10]. 

The control strategies mentioned previously used for the control of the DFIG requires the knowledge 
of the mathematical model of the system, the uncertainty due to the knowledge of the parameters of the system 
can degrade the performances of the system [11]. To solve this problem a robust and powerful control method 
called linear active disturbance rejection control (LADRC) is involved. This last one is based on the estimation 
and real-time compensation of internal and external disturbances, which can affect the system [12], [13]. 

The LADRC control strategy is based on the linear extended state observer (LESO) which 
constitutes its core and allows estimate and cancel the various external and internal disturbances of the 
system in real time [14]. This observer has the advantage that it does not require complete knowledge and 
exact system to be controlled. Furthermore, it allows to reduce the time to establish a detailed model of the 
studied system. 

In this present paper, the modeling strategy of DFIG is decribed; turbine, grid side converter (GSC), 
and filter, and then a robust control ADRC is applied to the control of the rotor side converter (RSC) and 
GSC. The test of the behavior of the system is highlighted using matlab/simulink software. Finally the last 
section concludes this work. 


2. RESEARCH METHOD 
2.1. Wind conversion system modeling 

The Figure 1 as shown in represent the different blocks of WCS. Actually, three-bladed horizontal 
axis wind turbine structures are the most widely used in the medium and large power wind turbine industry. 
These structures are based on the DFIG, and are equipped with two back-to-back type converters; RSC and 
GSC, the latter are controlled independently of one of the other. A capacitor is placed between the two 
converters to stabilize the direct voltage. This topology allows the converter to extract approximately 30% of 
the total power produced [14], [15]. 
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Figure 1. Block diagram of WCS based on DFIG 


2.1.1. Wind turbine model 
The function of the turbine is to convert the kinetic energy of the wind into mechanical energy of 
rotation. The extractable power available from the wind turbine (WT) is defined by (1). 


p.S. V3 


Paer = Cp. Pv = Cy. (1) 
Where V is wind speed, p is the air density and S = nR? is the area swept by the turbine. The power 
coefficient Cp expression is defined by (2) [13], [16]. 
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Where C, (A, B) is the power coefficient, A the speed ratio, and f the pitch angle. The gearbox ensures the 
connection between the turbine shaft to the generator shaft, its role is to adapt the speed between the turbine 
and the generator. The gearbox is modeled by a speed gain corresponding to the multiplication ratio given by 
the (4). 

dg 


G= (4) 


Qturbine 


In order to establish the evolution of the mechanical speed from the total torque applied to the rotor 
of the generator, the fundamental (5) of dynamics is applied [14]. 


düg _ T 
dt g 


J T Tem m f vg (5) 
Where f, is the viscous friction coefficient and J is total moment of inertia. The technique of extracting the 
maximum of the power consists in maintain the speed ratio A at its optimal value X opt which make the 
power coefficient to its maximum value cpmax [14]-[16], Figure 2 as shown in represent the block diagram 
of the turbine with maximum power point tracker (MPPT). 
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Figure 2. Block diagram of the turbine with MPPT 


2.1.2. Double fed induction generator model 
The dynamic model in a d — q reference frame of DFIG is writen by electrical (6) as shown in [15], [17]. 


ddsa 


Vsa = Rstsa + AE WsPsq 
he = Bain + 24 Oba 
Via = Rripa + S24 — oid, Mt 
Va = Roo TP + 0 Pra 
The electromagnetic torque Tom is illustrated by (7). 
Tem = 2D 2 (Psaira — Psaira) (7) 


To simplify the studied model, an orientation of the stator flux w along the d-axis is involved ($,4 = @,, 
sq = 0), also stator resistance R, is neglected [18], [19]. The DFIG stator voltages equations become (8). 


j dds 
Ms = Raise + a =0 (8) 


Voq = Oss 


Active and reactive stator powers are expressed by (9). 


Bulletin of Electr Eng & Inf, Vol. 11, No. 3, June 2022: 1210-1221 


Bulletin of Electr Eng & Inf ISSN: 2302-9285 C) 1213 


3 Vs Lm 
mes 2 Ls 74 
_ 3 Vss VsLm . (9) 
Que) T “TL, tra) 
Electromagnetic torque is given by (10). 
3 Lm š 
Tem = -3P T, sira (10) 


2.2. Double fed induction generator control 
2.2.1. Linear active disturbance rejection control model 

The LADRC control makes it possible to compensate for all real disturbances and model 
uncertainties of the system with the help of an LESO which is the heart of this control strategy [20], [21]. To 
illustrate the principle of the LADRC, a first order system with single input u and single-output y is 
considered as described by (11). 


VO = f(y,d,t) + bo- u(t) (11) 


dt ^ 
Where f (y, d, t) represents the dynamic of the model and the total disturbances (internal and external) of the 


system and bọ is the known part of the system. The representation in state variables of the system described 
by (11) is given as shown in (12) [20]-[22]. 


X,-f (12) 


In matrix form: 


(2)-G 96)* a+ Qe (13) 
y-a o(2) 


The first order LESO is illustrated in Figure 3, therefore the corresponding observer is presented by the (14). 


(B= Che 9) Gere " 
aa of 


Where ? is the estimated output, and L is the observer gain vector. For systems of n order, the general 
expression of the gains of the ESO observer are expressed by (15) [17]. 


_ (n+1)! i 
Boi = (n+1-i!i! ^? 


(15) 


Where i = 1,2,::+,n + 1 and n is the system ordr, and wọ denotes the observer's cut-off pulse. According to 
relation (15), the gains of the ESO for the Ist order system is expressed by (16). 


L=[Bo1 Boz F= [2«9 w ]* (16) 


It should be noted that a large pulsation makes it possible to have a good estimation of the states, 
however this risks increasing the sensitivity to noise. It can be determined as a function of the poles of the 
closed-loop system [23]. Figure 3 illustrates the block diagram of the first-order LESO observer. 

The estimated variables , — and £ĉ, f are used to compensate the effect of all disturbances by 
applying the following control law, given by the (17). 


uo-f 
bo 


(17) 
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Where: 


ug =k. (r — 9) (18) 


Figure 3. First order LESO 


Substituting (17) and (18) into (11) we get: 


uo-f 


by 


y= i + Bo: ff ter 9 (19) 


If f. ~ f and 9 ~ y we obtain: 

y =w =k(r—y) (20) 
Hence the following differential (21). 

_ytyer (21) 


We then obtain a closed-loop behavior of a first-order system having the pole: —k,. The controller 
gain —k, is chosen such that k, = w; = = where T, is the response time of the closed-loop system, w, is the 
closed loop natural frequency and r is the reference signal. The dynamics of the observer must be fast 
compared to that of the controller. Therefore, the poles of the observer are placed to the left of the poles of 
the closed-loop system. Generally, wo is chosen such that wọ = (3 ~ 10)w, [23], [24]. Figure 4 show the 
block diagram of the first-order LADRC. 
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Figure 4. Block diagram of first order LADRC 


2.2.2. Linear active disturbance rejection control for rotor side converter 

The simplified model of DFIG with stator field oriented shows that both direct and quadrature 
components of rotor current can be independently controlled, with each component having its own regulator. 
This control of the rotor currents makes it possible to generate the reference rotor voltages V4 ref and Vig ref 
to be applied to the pulse wide modulation (PWM) control stage of the converter [14]-[19]. The reference 
rotor current ig ref and tq per are given respectively by (22) and (23) [14]. 
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The rotor current ijj and i,q are controlled by LADRC controllers as shown in Figure 5. These 
expressions can be rewriten using the canonical form given by (24) and (25). 


dirg Ry. ` 1 
— = — — lra t Orlirg t+ — vV. 24 
dt oL, rd r'rq oL; rd ( ) 
dirg Ry . i Lm 1 

= — — lrg — WOrlra — Wr T — v 25 
dt oL Tq r'rd r oLsLy ds oLy rq ( ) 
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Figure 5. Control of the RSC using LADRC 


The LADRC controllers are then designed by (26) and (28) as shown in. 


dir 
2 = faC) + brota om 
Where: 
Rr, ; 1 
fal) = =z ira t aria, + Go — bro)Vra 
: ; (27) 


bo = a ug = V. 
ro oL ) ^d rd 
In addition: 


Ord = p C) + broug (28) 


dt — 


Where: 


= Rr , ` Lmós 1 
fac) z Tun "a g Wr (ira + 2r) + E 24 bro) Vrq 


1 (29) 
bro = 5 i Ug = Vr 


2.2.3. Linear active disturbance rejection control for grid side converter 

The control block of the GSC is shown in Figure 6, which includes a direct current (DC) bus voltage 
control loop to keep it constant and a current control loop in the filter (Rp, Lf) which is used to generate 
reference Veg yer and Vig yer voltages which will be applied to the input of the PWM command of the 
converter [25], [26]. 
The filtre currents expression can be writen using the canonical form given by (30) and (31). 


difa A, 
dt — 


Rf. ; 1 
E fd T Wslfq T Ly fa (30) 
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(31) 


Figure 6. Control of the GSC using LADRC 


Neglecting the losses in the filter (Rf, Lf), the active and reactive powers exchanged with the 
network are expressed by the (32) and (33). 


3 : 


Qr = AT (33) 


The power of the DC link capacitor C 1s described by (34). 


dUgc 
Pac m CUac "e (34) 
Where 
dUgc 1, ; 
a erc in (35) 
The power on the DC bus is given by the (36). 
Pac = P, — Pr (36) 


Where P., P; are the generator rotor side and filter side powers respectively. The power sign P. depends on 
how the generator works. P. is positive in hyper-synchronous mode (the rotor provides power to the grid), 
and negative in hypo-synchronous mode (the rotor receives power from the grid). The reference current 
iq ref comes from the internal voltage control loop of the DC bus voltage Uac. To impose a zero reactive 
power on the grid, the reference current i¢g_+e¢ can be given by (37). 


: 2 
lfa ref = zv, Qf ref (37) 
The canonical forms of LADRC are given by (38) and (40). 


di 
a = fraC) + broura (38) 


Where: 


_ Rf. ; 1 
fraC) = 1, f 7 Oslra t (= - bro) Vra 


f (39) 
bfo = ET $ Ufa = Vra 
In addition: 
diş = 
T = fral) + brotsa (40) 
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Where: 


_ 1 Rf., ; 1 
f = EM — irs — asia (7 bro) Va 


i (41) 
Po c1 4 Up Vg 


3. RESULTS AND DISCUSSION 

The WCS parameters used for the simulation are given in the Table 1 for the DFIG, Table 2 for the 
grid side and LADRC parameters. To justify the performance and robustness of the LADRC controller we 
tested and simulations of all parts of WCS, using the MATLAB/Simulink software. Two tests have been 
conducted under different conditions. In the first one the reference stator reactive power Qs ref has been set 
to zero and in the second one Qs ref has been set to —1e^ VAR. 


Table 1. The parameters of the DFIG 


Parameters Values 
Rated power Pn 1.5 MW 
Rated voltage Un 690 V 
Nominal frequency 50Hz 
Rated rotor speed 1750 tr/min 
Number of pole pairs 2 
Rotor resistance Rr 2.65 mQ 
Stator resistance Rs 2.63 mQ 


Stator leakageinductance Lso 0.1687 mH 
Rotor leakage inductance Lro — 0.1337 mH 


Table 2. Grid side and LADRC parametrers 


Parameters Values 
Rated power 1.5 MW 
Density of air 1.225 kg/m? 
Blade radius 30 
Gearbox ratio G 57 
Kp 400 
B01 2400 
B02 1440000 


3.1. Tracking test 

In this test, as seen in the Figure 7, a variable speed wind profile is applied to the WT. Figure 8 
illustrate the mechanical speed corresponding to the wind speed profile previously applied to the turbine. As 
seen in the Figure 9 the evolution of the power coefficient Cp(A,f) as a function of A and for different values 
of B. It can be observed that when the pitch angle p increases, the coefficient Cp decreases, and which results 
in a reduction in the wind kinetic energy captured by the turbine. The maximum of power coefficient Cp 1s 
given for the pitch f fixed at 0. The Figure 10 presents the evolution of the power coefficient in the case of 
the pitch 6 equal to 0. The aim is to extract maximum power and to obtain a maximum power coefficient Cp. 


212 MM 

En 3 E 1500 Y 
s" 5 o 1000 

a” 33 
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a " 1 2 3 4 5 a 1 2 3 4 5 

Time (seconds) Time (seconds) 
Figure 7. Wind speed profile Figure 8. Mechanical speed (rpm) 


The extracted mechanical power is shown in Figure 11. The Figures 12 and 13 illustrate respectively 
the evolution of the stator currents and rotor currents generated by the wind system. Figures 14 and 15 shows 
that the active power P; and the reactive power Q, at the stator level follow their references, proving that the 
control proposed by LADRC is robust and efficient. In Figure 16 the voltage Vdc of the DC bus is kept 


Proposed robust ADRC control of a DFIG used ... (Abdeslam Jabal Laafou) 


1218 O ISSN: 2302-9285 


constant with zero overshoot while Figure 17 represents the stator voltage Vsa and the stator current Isa. We 
can observe from Figure 17 that the voltage Vsa and the current Isa are in phase opposition which means that 


the DFIG operate in hyper synchronous mode. 


efficient Cp 


6 
Tip Speed Ratio 


Figure 9. Evolution of the Cp according to the speed ratio A and the pitch angle B 
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3.2. Comparison with backstepping controller 

To prove the robustness of the LADRC controller, a comparative study with the backstepping 
controller was established. Figures 18 and 19 represent the results of the active and reactive powers injected 
into the grid by the rotor. Therefore, it can be concluded that the power factor is equal to unit since the 
reactive power has been set to zero. From the previous results, it can be said that the LADRC controller is 
fast, robust, and stable than the backstepping controller. 
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Figure 18. Active rotor power exchanged with Figure 19. Reactive rotor power exchanged with 
the grid the grid 


3.3. Robustness test 

The principle of this test consists in changing the internal DFIG model parameters, to test robustess 
of the proposed control against the effect of the dynamic behavior of the system. Figures 20(a) and 20(b) 
show the evolution of the active and reactive power injected by stator side. The same wind speed profile was 
applied tothe turbine and the reference reactive power is chosen as shown in: 1) from t=0 to t=2.5 s and from 
t=4 s to t=5 s: Os-ref-0 VAR, ii) from t=2.5 to t=4 s: Qs-ref=-1 eSVAR. 

In this test, the rotor inductance Lr was changed by 2096 and the resistance Rr was increased by 
100%. One can clearly see in the Figure 20 that the variation of both resistance and inductance of the rotor 
has no effects when using LADRC control, which makes it possible to compensate the internal disturbances 
of the system. On the other hand and to demonstrate the performance of the LADRC regulator with respect to 
the backstepping controller, a comparison was made with a change of 35 96 on the rotor inductance Lr. 
Figures 21 and 22 show that the LADRC controller is more stable and robust against model uncertainties that 
affect system stability. 


25 3 35 4 45 5 Time (seconds) 
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Figure 20. Evolution of the (a) stator active power and (b) stator reactive power 
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Figure 21. Stator reactive power of the DFIG Figure 22. The electromagnetic torque of the 
for 35% Lr DFIG for 35% Lr 
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4. CONCLUSION 

In this paper, a linear control strategy based on LADRC controllers was established to control the 
RSC, GSC, and DC bus voltage in order to extract the maximum power regardless of the speed of the wind 
applied to the turbine. A judicious choice of the gain and the parameters of the LESO used in the LADRC 
regulation loop made it possible, according to the simulation results on the MATLAB/SIMULINK software, 
to compensate all modeling errors and attenuate the impact of varying system parameters. To prove its robust 
foundation, a comparison with backstepping control has been drawn. The results obtained showed that the 
LADRC controller is more robust due to its ability to estimate and compensate internal and external 
disturbances affecting the system under study. 
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